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Abstract 

We present a reanalysis of direct CP violation in the decay K L — > tt + tt~j. We point 

out an existing discrepancy between the theoretical and experimental definitions of e' H . 

Adopting the experimental definition of e' H , we estimate that e' H /e could be as large as 

a few times 10 both within the standard model and beyond. We discuss these estimates in 
detail and we also show how a judicious choice of E* cuts can increase the sensitivity of the 
observable e\ to the underlying CP violation. 
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1 Introduction 



The origin of CP violation remains an unsolved problem in particle physics. An important piece of 
this puzzle has been provided recently by the KTeV and NA48 collaborations that experimentally 
confirmed the presence of direct CP violation in the K — > tttt neutral-kaon decays. Interestingly, 
the measured results were larger than most existing standard model estimates at the time. The 
difference could be attributed to new physics, although updated analyses of theoretical uncertain- 
ties indicate that there is no serious disagreement with the standard model This situation 
underscores the importance of observing CP violation in other reactions. 

The decays K L s — > 7r + 7r~7 have long been recog nized |, |, |, |, [7j for their potential to test CP 
violation. The amplitude for each of these decays is conventionally divided into two contributions: 
inner bremsstrahlung (IB) and direct emission (DE). The IB is completely determined by the 
K — > 7r + 7r~ process that underlies it, whereas the DE part encodes additional dynamical features. 
The K s — ► 7r + 7r~7 decay is known to be dominated by the IB, but in the K h decay the IB is 
suppressed because the underlying K L — > 7r + 7r~ process is CP-violating. This makes the more 
interesting DE term in the K h decay more accessible, thereby raising the possibility of observing 
new direct CP violation in this mode. 

The presence of CP violation in K L — ► 7r + 7r~7 has been observed in recent experiments 0, |9|], 
and this is characterized by the parameter r)+- 7 , defined in analogy to the parameter rj + _ in the 
K — > 7r + 7r~ case. However, the result has not yet reached the level of sensitivity needed for detecting 
direct CP violation in K L — > 7r + 7r - 7. On the theoretical side, studies on direct CP violation in this 
decay have been performed by various authors ^ [7], [K| |Tl|, [12|] . In this paper, we revisit this subject 
for the following reasons. First, the existent experimental analysis of the CP-odd observable r?+_ 7 
is not consistent with the theoretical definitions because it assumes that this quantity is a constant 
over phase space. Second, KTeV will perform a new analysis of this mode with their recent data and 
it is therefore timely to update the theoretical expectations both within and beyond the standard 
model. Finally, we examine the possibility of measuring new, direct, CP violation in this mode 
from an analysis of the decay distribution JHJ without K s — K L interference. 

In the following section we introduce our notation for the relevant decay amplitudes and the 
CP- violating parameters. In particular, we introduce the parameter, e, to characterize the new 
direct CP violation in these modes. In Section |3|, we relate these parameters to CP- violating 
observables occurring in the interference between the amplitudes for K h and K s decays into 7r + 7r~7. 
By carefully treating the energy dependence in the amplitudes, we derive a relation between e, and 
the experimental observable e' + _ 7 . In this manner we extract the current limits on e. In Section 
we consider detecting direct CP violation in the interference between the inner bremsstrahlung and 
direct emission in the electric amplitude of K L — > 7r + 7r~7. Using a recent KTeV result, we extract 
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a bound on e. In Section |5] we present an estimate for e in the standard model. Finally, in Section [!], 
we consider potentially large new physics contributions to e that may arise in left-right symmetric 
models \TM and in generic supersymmetric models [ITH . 



2 Amplitudes and Direct CP- Violating Parameters 

The amplitudes for K° — > 7r + 7r~, 7r 7r° are conventionally written as 

M K0 ^ +n - = A e< + ±A 2 e< , M K o^ = A e< - V2A 2 e i5 % , (1) 

where Aj is the component for the tttt state with isospin J, and 5j is the strong-rescattering phase 
for the tttt state with angular momentum J. From K — > tttt data, one can extract in the isospin 
limit ReA ~ 2.72 x 10~ 7 GeV and u = Re A 2 /Re A ~ 1/22.2 |]|.[] The physical states K L and 
K s are given by 

\ k l) = / > l^s) = / > ( 2 ) 

/ 2 + 2e 2 ,/2 + 2e 2 



where e corresponds to CP violation in the kaon-mass matrix |~8| and we use the convention 
CP\K°) = — \K°). The CP-violating parameters in K — > tttt decays are 



ReA ' x/2^Re^ 2 ReA oy l ' l J 

corresponding to indirect and direct CP violation, respectively. It is also conventional to define the 
ratio of amplitudes 

r?+- = -r-^ = e + e' . (4) 

Experimentally it is found that 1 = (2.285 ± 0.019) x 10 -3 and that its phase is <p + _ = 

43.5°±0.6° ||16|| . The recent measurements of direct CP violation |l[ lead to a new world average [|l]| 
e'/e = (19.3 ±2.4) x 10" 4 . 

In the K — > 7T7T7 decay, the amplitude is generally decomposed into electric and magnetic terms. 
The electric part E receives contributions from both inner-bremsstrahlung and direct-emission pro- 
cesses, whereas the magnetic part M arises exclusively from the direct-emission. Our notation for 
the amplitude for K — > n + (p + ) n~(p_) 7(g), where K is any neutral kaon (K°, K°, K^, Ks), is 

4 (P+ +P-Yv - (p, -p_Yz 



m 



K 



+ M(K) , (5 ) 



m K 



1 Isospin violation significantly complicates this analysis 117 
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with the kinematic variables 

2k-(p + -p_) 2k -q 

v = 2 > z = — — ■ ( 6 ) 

m z K m z K 

In the kaon rest-frame, v = 2(E w+ — E n -)/m K and z = 2E*/m K , with E* being the usual notation 
for the photon energy in this frame. 

For K° — > 7r + 7r~7, the IB term, 

E m (K°) = A6Mk :-X , (7) 



7T + 7T . 



is completely determined [EH] by the amplitude for the underlying nonradiative decay K° 
Experimentally, it is separated by fitting the characteristic bremsstrahlung spectrum, which behaves 
as dT lB /dE* ~ 1/E* as E* — >• 0. Following Ref. 0, we write, guided by dimensional analysis, the 
DE terms 

E DE (K°) = elGy/^^z) , M(K°) = e \G 8 \ £ i£ M (v, z) , (8) 

so that the dimensionless form-factor £ E (^ M ) in the electric (magnetic) amplitude is expected to 
be of order one. Also, f n ~ 92.4 MeV is the pion-decay constant and G 8 is defined by 

A = V2GJAm 2 K -ml). (9) 

The corresponding amplitudes for K° — > 7r + 7r _ 7 are obtained by requiring CPT invariance, and 
they are 



E 1B (K C 



-4eM* K o_ 



m K (z 2 - u 2 ) ' ( 10 ) 

E DE (K°) = -e\G a \g&(-u,z), M(K°) = +e |G 8 | i&(-u, z) , 

where in A4* K ^ n7T and £g M the complex conjugation refers only to the weak CP- violating phases 
and not to the strong final-state interaction phases. 

The form factors £ EM can be expressed in a multipole expansion J3J. There is almost no exper- 
imental information on the electric form-factor, as it generates small corrections to the amplitudes 
in both K$ and decays. For this reason, it will be sufficient for us to assume that £ E is saturated 
by the leading multipole, El, and use 

&(!/,*) = F B e# , (11) 

where F E is a dimensionless complex constant, expected to be of order 1 by dimensional analysis. 
The final state interaction phase is 5{, reflecting the fact that in this amplitude the two pions are 
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in an I = J = 1 state. The magnetic form-factor, on the other hand, has been experimentally 



studied in some detail [^T]. £2] and found to depend mostly on z. For the remainder of this paper, 
we will take 



£mO>z) = £ M (z) e 



(12) 



corresponding to the leading multipole, Ml, with an E* dependent form-factor. In what follows we 
neglect any AI = 3/2 contribution to £em- 

For the physical kaon states K s and K L the IB amplitudes are given by 

AeM, 



E m (K s 



Kg— >7r+7T" 



m K \z* — v 



2\ ' 



E m (K L ) = V + -E m (K s ). 



The electric DE amplitudes for K h s — > 7i + n 7 are 



where 



& = V2ReF E 



e + i 



E 5 



Im F E Im A 
ReF E ~ ReA 



^ = V2 Re F E e i5 i . 



(13) 

(14) 

(15) 
(16) 



We have, as usual, dropped terms quadratic in weak phases. Similarly, for the magnetic amplitudes 
we have 



where 



M(K h>s ) = e\G a \£i$f(z) 



(17) 



SM 



e + i 



'Im£ 



M 



Im An 



Re £ M Re A c 



(19) 



It is then possible to define two CP- violating quantities associated with the electric and magnetic 
amplitudes. For the electric amplitude, one has the usual ratio || 00 



E m (K L ) + E DE (K L ) 
E 1B (K S ) + E DE (K S ) 



l+ - 4ReA ; 
i+ 4ReA V2 [ ] 



(20) 



2 Notice that we call it 77_| instead of ^ to differentiate it from the experimental observable. This is an 

important distinction that has not been appreciated in the literature. 
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Noting that 

\G&\ft™K _ L_ 

4A 4v/2 



m 



0.03 (21) 



K 



and that z 2 — v 2 is always less than 1, we write 



^ " ^ + 4ReA v/2 1 J 

^ |G 8 |/>^ReF E /We _ Im^ \ 2 _ ^ (22) 

/+ 4ReA V ReF E Re A o ; V ; V ; 



It is, therefore, convenient to define the parameter 



| G 8 1 f 2 m K Re F E ( Im F E Im A 



o 



4ReA \ ReF E ReA 0/ 



(23) 



to characterize new direct CP violation in K — > 7r + 7r 7. Thus, we have 

e^_ 7 = r? + _ 7 - = ee^ (^ 2 - v 2 ) , 5 e = b\ - 6° + vr/2 . (24) 

The quantity e' + _ 7 varies across the Dalitz plot, and for this reason we prefer e as a more natural 
measure of direct CP violation in this mode. 



Similarly, for the magnetic amplitudes, we can define 

M = M(K S ) _ M 



where in this case the direct-CP-violating parameter is 

and we have again dropped terms quadratic in weak phases. It appears that it is not possible to 
measure this parameter in an experiment that does not detect the photon polarization, as we will 
see in the next two sections. 

It was argued in the past J5], [TIJ that e' + _ 7 could be several times larger than e' because it 
predominantly arises from the interference between AI =1/2 components of the IB and DE terms 
and, therefore, it is not suppressed by the factor uj = Re A 2 /R,e A . Rather, since the IB and DE 
contributions are generated at orders p 2 and p 4 , respectively, in chiral perturbation theory, e' + _ 7 
has a suppression factor of p 2 /A 2 ~ mj^/A 2 , where A ~ 1 GeV. With the additional assumption 
that the weak phases in e' + _ 7 and e' are comparable, the enhancement was expected to bef\ e/e' ~ 

3 The original argument was done in terms of which is related to e by Eq. (p4|). The kinematic dependence 

of is such, however, that a considerable suppression results when it is integrated over phase space, as we will 

see in the next section. 
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{m 2 K /K 2 )/uo ~ 5. Using Eqs. (§ and (|23|), we find 

\G 8 \f'm K \F E \ f n \F E 



2V2u Re A 4c; m 



(27) 



A' 



From this result, one can see that if F E is of order one as we have speculated, then e ~ e' and 
there is no enhancement in this reaction. Equivalently, for the original estimate e ~ 5e' to be true, 
one would need F E ~ 5 . It is really impossible to distinguish between these two scenarios solely 
on the grounds of dimensional analysis. As we will see in Section [5], an estimate of F E based on 
factorization of the leading current- current operator in the effective weak Hamiltonian supplemented 
with vector- meson saturation of the p 4 strong counterterms suggests that F E ~ 1.7 and e is not 
much larger than e' within the standard model. 

A limit for ReF E can be obtained from the measured K s — ► 7r + 7r~7 rate, as a nonzero ReF E 
would give rise to a difference between the measured rate and the calculated IB rate [|3| . Including 
the DE electric-amplitude from Eqs. (|14"D and fll6|), and neglecting the magnetic contribution, we 
derive^ 

BR(K S -> tt+tt^) ~ 1.75 x 10~ 3 + 1.07 x 10~ 5 Re F E + 2.81 x 1(T 8 (Re F E ) 2 (28) 



for E* > 50MeV. Since measurements give BR(K S -> 7r + 7r~7, E* > 50MeV) = (1.76±0.06) x 



10 3 , we extract ReF E ~ 1 or —380 using the central value. Dropping the second solution, which 
is unnaturally large, we obtain 

Re F E = 1 ± 6 . (29) 

If the PDG's number @ BR(K S -> k + k--i,E* > 50MeV) = (1.78 ±0.05) x 10~ 3 is used instead, 
the result is similar, ReF E = 3 ± 5. Therefore, experimentally the question of the natural size of 
e/e' has not been resolved, and a value e/e' ~ 5 as in the original dimensional analysis estimate is 
possible. 



3 Interference between K L s —> tt + tt 7 amplitudes 

The parameter e may be measured in an experiment studying the interference between the ampli- 
tudes for K L and K s decaying into 7r + 7r~7. Such an experiment typically || § employs two K L 
beams, one of which is passed through a "regenerator", which coherently converts some K h to K s . 
It follows that the initial kaon state can be expressed as a coherent mixture 

\K h )+p\K s ), (30) 

4 In the first term (IB only) we have incorporated the complete K s — > 7r + 7r~ amplitude as extracted from 
data , and in the second term we have used an energy-dependent S® — S\ from Ref. [£4| . 
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up to a normalization constant, where p = \p\ e 1 ^ is the regeneration parameter. Then the number 
of decays into 7r + 7r~ 7 per unit proper time r is given by H 



dN 
dr 



{i/flX^e-As + r^_ 7 e-^ + 2Re [p 7 ^ s e iAmT ] e-^+^Jr/a} ( (31) 



where r L (r s ) is the (i^ s ) lifetime, Am is the K L -K S mass difference, is the partial 

width of K — > 7r + 7r~7, and 7 LS is an integral containing the interference between the K L and K s 
amplitudes. It is important to keep in mind that the quantities r^-^ and 7 LS are not integrated 
over all phase space. They depend implicitly on the cuts that define the region of phase space under 
study, and in the following discussion it is implied that all the terms in Eq. fl3~T|) are subject to the 
same set of kinematic cuts. 

We now proceed to examine the three terms in Eq. ([IT]) in detail. To this aim, we will make use 
of the following relations: 

1-^^+71—7 = \V+- | 2 ^K^it+tv-'y ) (32) 



which follows from Eq. ([13]); the ratio [22 



pM 

/ = ^ 0,85 , < 33 > 

which is determined experimentally from a fit to the decay spectrum that assumes that the inter- 



ference between the IB and an El-DE is negligible; and Eqs. (P0|), (pq). Consequently. 



pM 

r = r ^ +?r ~ 7 ~ 2.16 . (34) 

If the photon polarization is not observed, there is no interference between the electric and 
magnetic amplitudes. In this case the K s — ► 7r + 7r~ 7 rate can be decomposed into the sum of 
electric and magnetic rates 

pE _|_ pM 

Since the second term is CP- violating, it is convenient to rewrite it in terms of Eq. (|25| ) asQ 

p _ pE , I M |2pM 

Kg— >7T+7T 7 iCg— >7T+7r - 7 I'/H — 7I >7r + 7r_ 7 



;i + \v + -\ 2 Iv+^fr) r| s _ + ^ 7 ~ r| s _ + ,_ 7 . (35) 

The interference term in Eq. ( |3~1~1) can be written in the kaon rest-frame as 

7ls = J d[PS] {[E m (K L ) + E DE (K L )] [E; b (K s ) + E^(K S )} + M(K h ) M*(K S )} , (36) 



5 We have here assumed ?7+_ 7 to be a constant, which is appropriate for our purposes. 



where 



d[PS] = dcosfldE. 



(/3E;) 3 S in 2 g / _ 2^; 

327r 3 m^- \ m 



1 



(37) 



A' , 



with 6 being the angle between p + and q in the tttt rest-frame, and (3 = y 1 — 4m 2 / (m 2 ^ — 2E*m K ). 
We remark that d[PS] contains not only the phase-space factor, but also factors that resulted from 
summing over the photon polarizations and from contraction of tensor forms. Making use of the 
definitions in Eqs. Q21J), fl3J) and defining^ 



1 



A, 



/d[PS]eV 7 \E m (K s ) + E DE (K S ) 



(31 



one finds that 

7ls = K~ + 4- 7 + V+-,r [|r ?+ _| 2 + 2Re(r / ;_6 / + _ 7 )]}r i 

having dropped terms suppressed by additional powers of e or e. 
Turning now to the second term in Eq. (|31~D, we have 

"p 

A T — >7r+7T"~7 



(39) 



J d[PS] (|^b(^l) + ^de(^l)! 2 + |M(Jf L ) 



where the photon polarizations have been summed over. Then, from Eqs. ( p0| 
neglecting terms of orders e 2 and e 6 , we find 



(40) 
and (1331), 



'l + r) 



|^ + _| 2 + 2Re(r ? ;_e / _ 



Ac,^7r+7r~7 



(41) 



Collecting all these results, we can now rewrite Eq. ([TO) as 
dN 



dr 



oc 



+ 2 



77 + _| 2 + 2 Re (<_e' + _ 7 )] (1 + r) e"^ 



(42) 



|p| cos 



where 0^ is the phase of + e' + ). This rate equation can be used to extract e from measure- 
ments. The most recent experimental study || on these decays starts from the definition 



dN n s r K +7r _ 



dr 



\P\ 



{|p| 2 e- r/r s + |7/ + _ 7 | 2 (l + r)e- r ^ 
+ 2 |r? + _ 7 | |p| cos (Amr + p - e -(V- L +i/-s)/2} , 



(43) 



'Notice that both r Ks _ >7r+7r _ 7 in the denominator and the integral in the numerator of Eq. ( |38|) depend on the 



E* cut. 
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where N s is the number of K s regenerated. Comparing Eqs. and (f43|), we see that 



+~1 



(44) 



if terms of order (e/e) 2 are neglected. 

With a minimum photon energy, E* > 20MeV, we integrate over phase space ( — 1 < cos# < 1 
and 20MeV < E* < m K /2 - 2ml/m K , ) and use E{K S ) ~ E lB {K s ), to find 

e' ~ 0.041 ee i5 < , (45) 



where we have also used BR(iT s -> Tr+Tr^,^* > 20MeV) = 4.87 x 10~ 3 from Ref. |21[ and 
e 2 /(47r) = 1/137, and assumed 5 e to be constant. We have checked that a phase which varies 
with energy would not alter this result in a significant way. With an energy-dependent formula for 
5q — 5\ from Ref. f24|, we get e' + ~ (0.014 + 0.039 i) e, which is not different from Eq. 



we get e' + _ 7 — (0.014 + 0.039 i) e, which is not different from Eq. ( |45[ ) if 
5q — 5\ takes its average value of 17.4°. 

We see that the observable e' + _ 7 is suppressed by a factor of 0.041 with respect to e. This factor 
can be understood as the ratio of / d(phase space) / (z 2 — z/ 2 ) to / d(phase space) / [z 2 — v 2 ) 2 , which 
are the corresponding forms for the interference and IB terms in the rate. It is interesting to notice 
that the suppression factor decreases as the cut in E* is increased. This is due to the fact that 
with higher cuts the ratio of the IB contribution to that of the electric DE becomes smaller. For 
example, we find 



0.090 ee l4 « for E* > 50 MeV, 



(46) 



by means of BR(K S -> tt+tt^, E* > 50 MeV) = 1.76 x 10~ 3 from Ref. |2|. Clearly if the E* cut is 
increased much further, one has to reconsider the assumption of IB dominance in the K s amplitude. 
Experimentally, of course, an increased E* cut results in a smaller event sample. 



The best measurement (E* > 20 MeV) to date [7^ 



-+-7 



0.041 ±0.035 



(47) 



cxp 



translates into the one-sigma limit 



< 1.9 



Then, if we assume that there is no cancellation between the two phases in Eq. 
that lmF E could have a magnitude as large as 



(48) 

']), this implies 



Umi^ 



0.1 . 



(49) 
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The KTeV experiment is expected |^5| to reduce the uncertainties in |e' + _ 7 /e| to the 0.4% level, 
which would improve these bounds by a factor of 10. 

To end this section we comment on the observability of the quantity ?7+_ 7 for CP violation 
in the K s magnetic-amplitude. It is evident from Eqs. fl3"5p and ( |55| ) that ??+_ 7 occurs amongst 
terms of order e 3 or higher. Therefore, it does not seem possible to measure the direct-CP- violation 
parameter e+_ 7 , defined in Eq. (^6|), in this kind of experiment. 



4 Interference between IB and DE in electric amplitude 

The direct CP violating parameter e can also be detected in principle by analyzing in detail the 
decay distribution in K h — > 7r + 7r - 7 . Let us first consider the photon-energy (E*) spectrum used 
to separate IB from DE contributions. The E* distribution is given by 

dT K ^ n+w - y ((3E;) 3 



dEZ 32tt 3 



7 o*v. m K 



3 




l--—L\l sin 2 fldcos0 (\E lB (K L ) +E BE (K L )\ 2 + \M(K h )\ 2 ) . (50) 



Recent experimental studies |2T], |22[ of the E* spectrum assume that the DE process is purely 
magnetic and parameterize it with a form-factor-modified Ml contribution. By parameterizing 
the CP- violating El component in the DE electric amplitude as in Eq. (O) we can calculate its 
interference with the IB contribution. The presence of this term modifies the shape of the (E*) 
spectrum, and this generates a limit on e. 

We present in Fig. [I] an example of such a deviation. For this figure, the magnetic amplitude is 
parameterized in the form used by recent experiments [[H], [22 



^ " [m*-mC+2E*m K + ^ j ^ ' (51) 

where m p is the p- meson mass, and a l 2 are parameters obtained from the experimental fit. The 
solid curve describes the combination of the IB and magnetic DE processes, without any electric DE 
contribution. For this curve, we employ a 1 /a 2 = — 0.729 GeV 2 from Ref. |22j and adjust the value 
of a 2 in order to get a K L — > 7r + 7r~7 rate that matches the measured one (thus here a 2 ~ 3.09). 
The upper dashed-curve includes the IB, the magnetic DE with a 1 /a 2 as before and an IB-E1-DE 
interference using ReF E = and ImF E = 0.5 for illustration. In this case we have used a 2 ~ 2.40 
to keep the total rate fixed. The lower dashed-curve corresponds to the interference between the 
IB and El DE only. 

From Eq. (|5T)D, we find the interference between the IB and DE components of the electric 
amplitude to be 

r^ +7r - 7 = / dpS] 2Re[E: B (K L )E DE (K L )] , (52) 
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40 



80 120 

E* (MeV) 



160 



200 



Figure 1: Photon-energy (E*) distribution in K L — ► 7r + 7r~7. The solid curve represents the sum of 
the inner bremsstrahlung and the magnetic direct-emission contributions. The dotted curve is the 
IB alone. The upper dashed-curve corresponds to IB plus Ml-DE contributions as in the solid curve, 
plus an IB-E1-DE interference with F E = 0.5i, as explained in the text. The lower dashed-curve 
shows the interference between the IB and El DE only. 



with the integration ranges being —1 < cos 9 < 1 and 20 MeV < E* < m K /2 — 2m\/m K . Using 



Tk^^-^E; > 20 MeV) = 5.86 x 10" 19 MeV 
from Ref. [ 16*f , Eqs. (^) and flTBD , as well as 5q — 8\ = 17.4° in 5 e , we obtain 



(53) 



iE,int 



5.0 x 10 3 |r/ + _|ecos(5 e -0 + _) + 1.8 x 10 2 \er) + _\ ReF E cos(5 e - tt/2) 



0.023 



19 x 10~ 4 ReF E . 



(54) 



An analysis of KTeV data with the same E* cut, gives [p~3 1 



-iE,int 

Kj — >7T+7T~ 7 



< 0.30 (90% C.L.) . 



(55) 



Taken literally, this bound implies that, 



+ 0.04 Re F F 



< 13 



(56) 
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and assuming that the two terms do not cancel, 

|ImF E | < 1, | Re F E | < 350. (57) 

These limits are much weaker than Eqs. (^8|), (|49"D, and (p9|). However, it is important to notice 
that the bound, Eq. (|55"1), is not particularly strong. In fact, it is amazing that this interference 
between the IB and El-DE contributions, which is expected to be small, could account for as much 
as one third of the experimental rate. This is due in part to the fact that an analysis of the full E* 
distribution is not the optimal observable to isolate this term as we can see in Fig. [l]. 

It may be possible to improve the bound, Eq. (^) by restricting the analysis to the region of the 
E* distribution where this term is most important. A glance at the dashed curves in Fig. p] suggests 
that the cuts 50 MeV < E* < 90 MeV, for example, would improve the bound. With these cuts, 
we find 

pE,int 

x L ^ +7r - 7 ^ 0.029 1^7 + 1.1 x 10~ 3 ReF E , (58) 

where we used 5$ — 8\ = 21.3° (the average phase difference in this region) and 

r^ L _^- 7 (50MeV < E* < 90 MeV) ~ 1.84 x 10~ 19 MeV , (59) 

obtained by including only the IB and magnetic contributions (since these should dominate the 
rate).[] These numbers suggest that this simple set of cuts could improve the bound on e by about 
30%. However, this procedure depends on the assumed E* dependence of Eq. ([51]), and of F E 
(independent of E*). 

A more detailed analysis of the decay distribution is desirable in order to separate all the 
contributions. Unfortunately, we find that the decay distributions for the Ml term and for the 
IB-El interference term are remarkably similar when the photon polarization is not observed. We 
illustrate this by presenting a cos 9 distribution in Fig. ^. We can see in this figure that the shapes 
of the Ml distribution (solid line), and IB-El interference term (dashed line) are very similar. The 
only significant difference between the two is that the interference term could be negative, depending 
on the phase of F E . 

5 Standard-model contribution 

In this section we review the estimate of F E that exists in the literature and we apply it to an 
estimate of e within the standard model. 

7 For the magnetic part, a l /a 2 — — 0.729 GeV 2 and a 2 — 3.09 have been used. 
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-1.0 -0.5 0.5 1.0 

cos 

Figure 2: Distributions of various components of the K L — > 7r + 7r~7 rate in terms of cos8, where 9 
is the angle between the ir + and 7 momenta in the irn rest-frame. The solid curve represents only 
the magnetic contribution. The dotted curve is the IB alone. The dashed-curve corresponds to the 
IB-E1-DE interference with F E = 0.5i. 

The SM contribution to the DE electric-amplitude in K L — > 7r + 7r~7 is generated at short 
distance by both four-quark and two-quark (s — > d r y ) dg) operators with known coefficients. The 
matrix elements of these operators arise at 0(p A ) in chiral perturbation theory. They are presumably 
dominated by the chiral realizations of the four-quark operators that transform as (8 L , 1 R ) under 
chiral rotations, as well as by 0(p A ) loop diagrams []7|, |26[. The contribution from the s — > dj,dg 
operators in the SM starts at order p 6 , and since the coefficients of these two-quark operators are 
quite small, it is appropriate to neglect them in the standard model estimate. 

The 0(p A ) weak chiral Lagrangian transforming as (8 L , 1 R ) is given by p 

4? = GJl £ NW + H.c. , (60) 

i=14 

where N t are dimensionless coupling constants and 

W u = i(A {F% + rfFpJ , IfV}) , W 15 = i(AL" (F%, + U^F%U)L V ) , 

(61) 

W 16 = i( A { F l - UtF > > L " LU }) > W 17 = i(AL" (F£ - tfF%U)L v ) . 

In these formulas, (•■■) = Tr(- ■ •), 

A = ±(A 6 - iA 7 ) , = F* = -eQF^ , U = e 1 ^ , = iC/+ D^U , (62) 
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where F = d ' A u — d v A„ is the photon field-strength tensor, Q = diag(2, — 1, — 1)/3 is the 
quark-charge matrix, <p is a 3x3 matrix containing the octet of pseudo-Goldstone bosons, and 
D U = d U + ie[Q, U]A . Under chiral rotations, these fields transform as 

F% - V h F%,Vl , F% - V*F%Vl . U -> V K UV£ , L M - • (63) 

In terms of this effective Lagrangian, the SM contribution to F E can be written as, 

Fe = V2P {NlA ~ Nl5 ~ Nw ~ Nir) + ° ip4) l0 ° P t6rmS ' (64) 

This combination of constants (iVj is not known from experiment so one needs to resort to model 
calculations as in Ref. f28| . The loop terms, on the other hand, are known [f2£fl . 

The standard assumption is that the real parts of the N i (i — 14, • • • , 17) arise mostly from the 
dominant octet quark-operator Q l — Q 2 (in the notation of Ref. f3H||). They have been calculated 
in a couple of models in Ref. with the result, 



Re (N u - N 15 - N ie - N 17 ) = ^ k f , (65) 

where kt — 1 corresponds to the naive factorization model and kf = 1/2 corresponds to the so-called 
"weak deformation model". Taking kf = 1 for illustration, and including the loop terms calculated 
by D'Ambrosio and Isidori (the 0.9 term), one obtains the estimate, 

|ReF E | SM = f - + 0.9 = 1.7. (66) 

2 V2 f^ml 

In order to estimate e one also needs an estimate for the imaginary parts of the matrix elements, 
or equivalently, of the N i . Within the standard model these imaginary parts are expected to be 
dominated by the gluonic penguin operator. Unfortunately a detailed estimate of the bosonization 
of the penguin operator at order p A does not exist, and we have to resort to dimensional analysis 
arguments which indicate that ImF E /ReF E is of the same order as Ini74 /Rev4 ~ \^2e'/u ~ 
1 x 10~ 4 . Taking ReF E = 1, and assuming that there is no large cancellation between the two 
terms in Eq. we arrive at the estimate 



< — T J^ I — -r ~ 3 x 10" 3 . (67) 

SM 2V2\e\m K ReA 



This number is almost three orders of magnitude below the current experimental limit in Eq. 
It should be clear, however, that this is no more than an order of magnitude estimate. 
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6 New physics and short-distance s —> dj, dg transitions 

Physics beyond the standard model modifies the short distance coefficients of the two and four-quark 
operators that contribute to the DE electric-amplitude in K h — > 7r + 7r~7, and may also generate 
further operators. In this section we discuss the s — > dj, dg operators, as they can be significantly 
enhanced in certain models [12|, [31], |32], p3| . 

The effective Hamiltonian responsible for the short-distance s — > dj, dg transitions can be 
written, following the notation of Ref . |3l| , as 



H eS = C+Q+ + C-Q- + C+Q+ + C-Q- + H.c. , (6* 



where are the Wilson coefficients and 



Qy = ^{sL%„d R ±s R a^djF^ , 

16tt 2 V J (69) 

Qt = ^(s L ^^ad R ±s R a^t a d L )G^ 



± 

1,9 



are the so-called electromagnetic- and chromomagnetic-dipole operators, respectively, with eQ a 
— e/3 being the <i-quark charge and G^ u being the gluon field-strength tensor. We notice that Q 
transform as (3 L ,3 R ) ± (3 L ,3 R ) under SU(3) L xSU(3) R transformations, and Q+ g (QZ g ) are even 
(odd) under parity. We further notice that (QZ g ) are also even (odd) under a CPS trans- 
formation (a CP operation followed by interchanging the s and d quarks). These magnetic-dipole 
operators are expected to produce the most important contributions to the s — > d'y, dg transitions. 
In this section we estimate the impact of these operators by employing a chiral-Lagrangian ap- 



proach combined with naive dimensional-analysis |34|, |35| to evaluate the necessary hadronic matrix 
elements. 

The chiral Lagrangian that represents the s — > d'y, dg operators is constructed as follows. We 
begin by defining the objects 

A LR = A RL = A , (70) 
with which we can rewrite in the form 

05 = ^5 (& A LR<V & ± fe V<V Ql) , (71) 
where q = (u d s) T . If we imagine that under a chiral rotation 

Alr - V L X LR V^ , A RL -> V R X RL Vl , (72) 

where Vl,r £ SU(3) L r , then the quark operators inside the brackets in Eq. ( [71] ) would be chirally 
invariant. Consequently, we construct the effective Lagrangian by writing chiral invariant terms 
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with one power of A LR or A RL assumed to transform as in Eq. (0), and with an explicit photon 
field strength tensor F^ u . We also require the effective Lagrangian to have the same parity and 
CPS structure of Q^. In this regard, A LR and A RL interchange places under parity, but remain 
unchanged under charge conjugation. 



There are many possible chiral realizations of Q*. As relevant examples at leading order, p , we 
write down 



ifieQ d (\ LK ULL±X RL LLU^)F^ + H.c 



(73) 



where have the same symmetry properties as those of Q^. Employing naive dimensional- 
analysis ]34], we obtain the order-of- magnitude estimate 



^7 fir fir 



P ± 

^ 16vT 2 A 



(74) 



where A = Aixf^. Notice that £w are not suppressed by light-quark masses, as is appropriate 
for the new-physics interactions of interest in Eq. (|68|). As was noted before, within the SM the 
short- distance operator is suppressed by light-quark masses, and this results in a different chiral 
Lagrangian (of order p 6 at least) involving the usual chiral-symmetry breaking factor |7j]. 

The contribution of the s — *> d'y operators to the direct-emission electric-amplitude of K° — > 
7r + 7r - 7 comes from C^~' and is given by 



\/2 f3_ m\ vr? K 



c: 



(75) 



Then, if the new-physics contribution to F E is such that ImF E /ReF E 3> lmA /ReA (i.e. CP 
violation is dominated by new physics) in e, the contribution of C~ to e/e is 





ImC- 




384^2 7T 3 





6.9 x 10 5 GeV 



ImC" 



This and Eq. 



imply the current limit 



ImC: 



< 2.7 x lCT 6 GeV" 1 . 



(76) 



(77) 



Future measurements [^| may improve this value by about a factor of 10. 

Turning now to the chiral realization of the s dg operators Qg, we rewrite them as 



Q 



± 



167T 



(7* 



Since do not contain the photon field, in constructing the corresponding chiral Lagrangian 
involving an F , we employ the chiral field-strength tensors Fh v and F^ u , which were defined in 
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Eq. flB2|). These tensors transform under parity as F^ R — > (F Rh ) and under charge conjugation 

Thus, examples of chiral Lagrangians at order p 4 that correspond to are 
4 ±} = ^ < A lr^ m + ^*^) V ± A RL L M (j£ + rtF^u) L v Tfl) + H.c. , (79) 
where, from naive dimensional- analysis @, 35] , 



0t 



167T 2 A 



^9 S 



10) 



For numerical estimates we use g s ~ v47r, corresponding to a strong coupling a s ~ 1. The 
contribution to F E arises from C g ~^ and is given by 



3 

(F ) 9s m K q- 



^1) 



We again assume that CP violation is dominated by the new physics, in such a way that ImF E /ReF E 
Imy4 /ReA , to obtain 



9s m K 


ImC; 




384V/2 7T 3 


e\f*ReA 



2.4 x 10 6 GeV 



ImC: 



52) 



which, with Eq. (f48"|), implies the limit 



ImC 



< 7.8 x 1(T 7 GeV" 1 . 



13) 



Since it is known |37]] that Q also contributes to e', we can use the measured value 



of e' to derive another limit on the contribution of Q g to e. The leading chiral realization of Q g 
that contributes to the K L — > txtx amplitude is of order p 2 , as the 0{jP) realization does not 



contribute once tadpole diagrams have been properly taken into account []36|]. An example that we 
can construct isPI 



4 2) 



%fH\(U-tf))(L»L\ + H. 



54) 



where 7 9 is a dimensionless coupling constant. Using the naive dimensional- analysis prescribed in 
Ref. f3"5j|, we find the order-of- magnitude estimate 



7a 



Cg Ids 
167T 2 



15) 



i A similar Lagrangian has been given in Ref. 



17 



The resulting amplitude is 



(A) 



275^ " 2m - 



(86) 



leading to 



10 



(m 2 K — 2m^ 



ImC: 



1.4 x 10 fc 



ImC; 



|e| ReA An 2 

Assuming that the current value \e'/e\ ~2x 10~ 3 is saturated by the Q~ contribution yields 



*7) 



ImC; 



< 1.6 x 10" 



(88) 



which is a much better constraint than Eq. fl33|). Furthermore, combining Eqs. ( |8T| ) and fl32| ) results 
in 



A' 



96^2 Trwe^ 2 (m 2 K -2ml) 



e' 




e' 




~ 1.8 




e 


9 


e 



< 4 x 10 



-3 



59) 



This is twice as large as the standard model result, but the two should be considered equivalent 
within the uncertainties of our estimates. From this we conclude that improved measurements of e 
are more important to place bounds on Q~ . Bounds on Q~ from e are not likely to be competitive 
with bounds from e'/e in the foreseeable future. We now turn our attention to two specific models 
for Q~ . 



6.1 Left-right symmetric models 



In left-right symmetric models the coefficients of Q~ and Q~ can be enhanced because the mixing 
of left- and right-handed IV-bosons removes the helicity suppression present in the standard model. 
Variations of this model have been studied in the context of b — > 37 in detail |39], flOfl . 

We start from the effective Lagrangian that results from integrating out the heavy right-handed 
W. This can be written down directly by following the formalism of Ref. [0. In the unitary gauge, 



r - tl2 



92_ 

V2 



U R C R 



WJ + h.c. 



(90) 



where V is a 3 x 3 unitary matrix having elements V , 



V qq ,K^ q ,, with V qq , being CKM- matrix 



elements and K qq , complex numbers. In writing C eS above, we have ignored modifications to the 
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left-handed I^-couplings which do not lead to enhanced effects. Using the results of Refs. [[39|, f|2 
we write, 



^~^F ^ — ^ * / R,* R, \ "^R,h(^o) 

Z^VqdVqa { K qs ~ K qd) m q 

V 2 q=c j 



(91) 



where x = m 2 Jm^ and 



RH 



-3x 2 + 2x , 5x 2 - 31z + 20 
lnx — 



fa; - l) 3 



Crh( x ) 



6x lnx 



3 + 3x 



6(x-l)2 • (a; -l) 3 (x-1) 2 (92) 

For our numerical estimates, we will use a s (m z ) = 0.119, m c = 1.25 GeV, m t = 173.8 GeV, and 
the CKM-matrix elements in the Wolfenstein parameterization from Ref. [44|: A = 0.22, A = 0.82, 
p = 0.16, and r\ = 0.38. This gives, 



°7,rh — 



i x io 3 v cd v; s («£ - «*) + 7 x io 4 v td v* 



K 



td 



x 10- ' GeV 



-i 



(93) 



c: 



It then follows that 



-4 x 10 2 V cd V c * s (4* - «£) - 2 x 10 4 (4* - 4)] x IO" 7 GeV" 1 . (94) 



7 ,RH 



15 Im (K**-K*)+Im [(1.3-0.6i)( 



_R* 



(95) 



These same couplings contribute to e'/e in this model through 



g,RH 



12 Im 



Im 



(0.8-0.41) («2T-k 



(96) 



From these results, we see that the contribution of to e is also constrained by the contribution 
of Q~ to e'. For example, if all Im(/c?f — K^ d ) are of the same order, then the c-quark contributions 
dominate both observables and |e| ~ 1.3 |e'| . Similarly, if the t-quark intermediate-state dominates, 
|e| ~ l-7|e'| . If one includes perturbative-QCD running of the Wilson coefficients from the W- 
mass scale down to a scale near the charm-quark mass []30], |43| , these numbers slightly change to 



|e| ~ 1.5 |e'| and 1.9 |e'|, respectively. With |e'/e| Sj RH ~ 2 x 10 3 , this results in 

~ 4 x 10~ 3 , 



(97) 



7,RH 



which is comparable to the contribution of Q g as in Eq. fl8"9p. In full generality, e/e and e'/e are 
proportional to different combinations of the k's, so that it is possible for e to be significantly larger 
than e', although this does not seem likely. 
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6.2 Supersymmetric Models 



In certain supersymmetric models, one can generate the s — > drj operators at one-loop via inter- 
mediate squarks and gluinos resulting in large C 7ff . The enhancement is due both to the strong 
coupling constant and to the removal of chirality suppression present in the standard model. We 
follow Ref. |TJ] and work in the so-called mass-insertion approximation. The full expressions for 
C~ susy can be found in Ref. |15[]. Here we are interested only in the terms enhanced by m g /m s , 
with m- being the average gluino-mass, and they are 

7ra B (m § ) 



^7,susy( m tf) 



2l)LR 



«)rl 



F (x ) 

susjK^gqJ 



(98) 



2l)LR 



2l)RL 



G (x ) 

K ' susy v^gqj 



X 



91 



where the <5's are the parameters of the mass- insert ion formalism, 
average squark-mass, and 

Ax (1 + Ax — 5x 2 + Ax lnx + 2x 2 Inx) 



m^/m^, with rriq being the 



susy 



x) 



3(1 



x 



^susyl 3 ') 



x (22 — 20x — 2x 2 — x 2 In x + 16x In x + 9 In x) 



(99) 



3(1 



x\ 



For our estimates, it is sufficient to approximate _F E 



susy 



X 



F (V 

susyV / 



^susyl 1 



2/9 and G. 

5/18 . This approximation introduces an error smaller than a factor of two. 
Then, by means of Eqs. ( |76D and (|87|), we obtain 



susy 



X 



7, susy 



92 



0.4 



500 GeV 



v a s (500 GeV) 
e 

9,susy 



m g 



Im 



d \ 

12JLR 



d \ 

12JRL 



(100) 



We find that the Q~ contribution to e' constrains the Q~ contribution to e as it also happened in 
the left-right model. Once again with \e'/e\ 
reach 

e 



9,sus y ~ 2 x 10 3 , the contribution of C 7susy to e can 



8 x 10 



-4 



(101) 



7, susy 



If we incorporate the running of the coefficients down to a low-energy scale |j3*I] , we will instead 
have |e/e 



7, susy 



91). 



1.7 x 10 , which is roughly comparable to the contribution of Q g in Eq. 

1 

The factor of three can be traced back to the use in Ref. 11211 of a hadronic matrix element 



Our estimate, translated into e' + _ 7 , is three times smaller than the estimate of this parameter in 

Ref. L _ 
three times larger than ours. These differences are a good indication of the level of precision that 
can be expected from this type of estimates. 
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7 Summary and Conclusions 



We have reanalyzed direct CP violation in the electric amplitude for the decay K L — > 7r + 7r~7. The 
dominant CP violating observable arises from an interference between the IB (inner bremsstrahlung) 
and the El direct emission amplitudes. The previous theoretical definition of ^+_ 7 results in a 
quantity that varies with kinematic variables and that differs from the previously used experimental 
definition of 77+- T . To clarify this situation we have introduced a new theoretical quantity e, to 
parameterize new direct CP violation in this decay. This quantity is a constant, and by simple 
dimensional analysis is expected to be between a few and five times larger than e'. 

The experimental observable e' + _ 7 is related to e by a normalized integration over phase space. 
We find that the very different E* dependence of the IB and El-DE amplitudes introduces a large 
kinematic suppression into e' + _ 7 - In particular, for E* > 20 MeV, |e' + _ 7 | = 0.041 |e|, and thus much 
smaller than e'. This kinematic suppression can be reduced by increasing the minimum E* accepted 
at the cost of diminished statistics. 

From presently available data we have extracted the bound |e/e| < 1.9, significantly larger than 
the theoretical expectation |e/e| < 0.01. 

We have estimated the value of e in several models. The estimates are hindered by unknown 
hadronic matrix elements and must be considered order of magnitude estimates. The general 
conclusions from these estimates are 

• If the CP- violating phase of the El-DE amplitude is similar in size to the CP violating phase 
of Aq, and the two do not cancel, then e is about as large as e' and could be as much as five 
times larger. This is a refined form of the naive dimensional analysis estimate. This is the 
case in the standard model estimate of e where the CP-violating phase cannot be computed 
and is simply assumed to be of the same order as the phase of A . 

• Beyond the standard model, in models where the chromomagnetic dipole operator is enhanced, 
the maximum value of e is directly constrained by the fraction of e' that is attributed to the 
new physics. Without invoking a fine-tuned cancellation in e', this implies that in these models 
e is also at most a few to five times larger than e'. This happens, of course, because the bound 
from e' is equivalent to the simple assumption that the CP violating phases of the El-DE 
amplitude and A are of similar size. 

• In principle, e could be larger in models in which it is induced primarily by an electromagnetic 
dipole operator which does not contribute significantly to e'. In the specific models that we 
studied, however, the coefficients of the electromagnetic and chromomagnetic dipole operators 
are highly correlated. For this reason, in these models, the maximum e is, once again, limited 
by the fraction of e' that we attribute to the new physics. 
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To conclude, e is expected to be comparable to e' although it can be as much as five times 
larger both in the standard model and beyond. Unfortunately, the experimental observable e' + is 
kinematically suppressed, calculated to be ~ 0.04 — 0.09e for E* > 20 — 50 MeV. 
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